Dipolar Poisson-Boltzmann Equation: Ions and Dipoles Close to Charged Surfaces 
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We present an extension to the Poisson-Boltzmann model where the dipolar features of solvent 
molecules are taken explicitly into account. The formulation is derived at mean-field level and can be 
extended to any order in a systematic expansion. It is applied to a two-plate system with oppositely 
charged surfaces. The ion distribution and profiles in the dipolar order parameter are calculated 
and can result in a large correction to the inter-plate pressure. 
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Charged objects (ions, interfaces and particles) im- 
mersed in liquids play a central role in electrochemistry, 
colloidal science and biology ranging from electrolyte ap- 
plications, stabilization of colloidal suspensions, protein 
folding and its biological activity, and even in protein 
aggregation [Elli, II g. 

The most commonly used model - - the Poisson- 
Boltzmann model (PB) [l|,H3| — assumes point-like ions 
immersed in a continuum dielectric media and treats the 
system in a mean-field approximation. The medium is 
modeled by a homogeneous and isotropic dielectric con- 
stant. This model is simple, elegant and efficient. It is 
in good agreement with experiments for monovalent ions 
up to energies of order of kgT. However, careful mea- 
surements of the forces between two charged surfaces at 
nanometric scale show strong deviation from the simple 
PB picture @. In particular, the assumption that the 
continuum dielectric medium is homogeneous does not 
take into account the strong dielectric response of water 
molecules around charges. The discrete moments of wa- 
ter molecules will orient themselves close to charged ions 
and surfaces giving rise to hydration shells and to hy- 
drophobic interactions, which can be measured at short 
distances, for example, between two charged plates (sur- 
face force balance apparatus) . These hydration phenom- 
ena are very important in many biological processes such 
as protein folding, protein crystallization and interactions 
between charged biopolymers inside the cell. 

Most studies other than the PB rely on one of sev- 
eral theoretical techniques. Monte Carlo (MC) @ or 
Molecular Dynamic (MD) Q computer simulations take 
into account the discrete nature of the dipolar molecules. 
A second approach relies on liquid state theory, integral 
equation and other methods jlOl . [Tlj . In simple planar ge- 
ometry the latter gives good agreement with the MC and 
MD simulations. However, all these methods are rather 
cumbersome and involve heavy computation resources. 
In addition, they lack the simple physical picture pro- 
vided by a Poisson-Boltzmann type of approach. 

In this Letter we propose another approach called the 
Dipolar Poisson-Boltzmann (DPB). Unlike the PB model 
where the solution is characterized by a homogeneous 
dielectric constant, in the DPB model we coarse grain 



the interaction of individual ions and dipoles interact- 
ing together. This makes the DPB an analytic extension 
of the PB formalism. Although it is done on a mean- 
field level, it includes some aspects of the discrete nature 
of the dipolar solvent molecules and how they modify 
the ion-solvent interactions. We show that such correc- 
tions to the PB treatment are important in predicting 
dipolar profiles close to charged surfaces and result in a 
strong deviation from their average value. Furthermore, 
the DPB model can, in principle, be expanded to any 
desired higher order in a systematic expansion. 

Consider a system composed of Nd mobile dipoles each 
with a dipolar moment p and / species of ions immersed 
in a continuum dielectric medium with a weak dielectric 
response (the justification for this system set-up is elab- 
orated below), e > £oj £o being the vacuum permittivity. 
Each ionic species has Nj ions of charge qje, j — I ... I, 
where e is the electron charge. In addition, the system 
includes a fixed charge distribution p/(r). The charge 
density created by a point dipole p at point ro is given 
by p,i(r) = — p- V5(r — ro). Thus, the total charge density 
is 

N d I Nj 

p(r) = -^p i .V ( 5(r-r l ) + ^^ ( z J e ( 5(r-R^) + p / (r) 

i—1 j—1 i—1 

(i) 

where denotes the position of dipoles of moment p^ and 
Rp- 1 are the positions of ions of type j. The canonical 
partition function is given by 

~~ N d \Il!_ N-l J i=1 PiUi^W" K * 

j — 1 J 
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where v c (r) denotes the Coulomb potential. Using a stan- 
dard Hubbard-Stratanovich transformation, 

Z = J Xty(r) exp (-^-J d 3 r [V0(r)] 2 
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where e = e^e r is the medium dielectric constant (in SI 
units) and (3 = l/T is the inverse temperature (where the 
Boltzmann constant, ks, is set to unity). The fugacities 
of the dipoles and i th ion species, A^ and Aj, respectively, 
are derived from the relations: Nd = A^gf-logZ and 

N i = X i ^\ogZ. 

Assuming that each molecular dipole has a fixed mag- 
nitude, |p| = po we sum now over the {p} degrees of 
freedom and obtain 
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The DPB equation is obtained as the saddle-point of 
the action (U|) [where we have used W(r) = i<fi(r) to denote 
the physical electrostatic potential] 
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and the function Q(u) — coshit/u — sinhit/w 2 is related 
to the Langevin function C{u) — cothw — 1/uhy Q = 
{smhu/u)C. One recognizes in ([5]) the usual terms of 
the Poisson-Boltzmann equation (the first two terms on 
the RHS), while the last term is the divergence of the 
polarization contributing to the induced charge density. 
The local polarization density (square bracket) in Eq. ([5]) 
is the product of the dipole density, sinhu/u, and the 
average dipole moment given (on a mean-field level) by 
the Langevin function. 

In the following we study a dipolar solvent with 1:1 
salt confined between two oppositely charged planes [l3j . 
While the dipolar effects are small for two similarly 
charged plates, they are pronounced for anti-symmetric 
plates and yield a spatial variation of the dielectric con- 
stant. Choosing the charge density to be =F<7 for the two 
plates located at z = id/ 2, the potential, ionic profiles 
and dipole density depend only on the z coordinate per- 
pendicular to the planes and ([5]) becomes 



e5l"{z) = -2c s esinh/3e# + aS(z + d/2) 

aS(z-d/2) + c dPo ^-\g(fip Q V) 
dz L 



(6) 



where we assume that the system is in contact with a 
reservoir containing a dipolar fluid of concentration Cd 
and salt of concentration c s so that Xd = Cd and X s = c s . 

The boundary condition at the z = —d/2 charged 
plane is 



= c d poG{(3poK)+cr 



(7) 



and the electric field E = —ty' is the same, for the anti- 
symmetric system, as on the other plane. Note that the 
usual Neumann boundary conditions for the PB equation 
includes now the polarization induced surface charges. 
We find that for strong enough surface charge densities 
the induced charge can be substantial and corresponds to 
a large modification of the standard boundary condition. 

From ^ we obtain the first integral which is equiva- 
lent to the contact theorem expression for the pressure 
difference II = P in — P out 
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This equation allows to express ^(z) as a of function 
fy' and thus solves © by a simple quadrature. The first 
two terms in II are the usual PB contributions, the first 
being the electric field and the second the mixing entropy 
of the ions. The other two terms are the specific terms of 
the DPB model. The first is the enthalpic contribution 
related to the orientation of the dipoles in a local elec- 
tric field. The last term is the rotational entropy of the 
dipoles. The pressure at any point z is calculated with 
respect to the pressure exerted by the bulk reservoir out- 
side the plates. 

Another way to interpret is to write it as a PB 
equation with an effective field-dependent dielectric con- 
stant e cS (E) = e ef(E) replacing the e on the LHS. The 
non-linear dielectric response is given by 

e eS (E)=e+^g(0p o E) (9) 

For weak fields one can expand the function Q to first or- 
der and obtain the standard PB equation e cff, I'"(z) w 
2c s esinh/?e x I' with an effective homogeneous dielectric 
constant 



,eff 
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(10) 



This result for dielectric response of molecules with in- 
trinsic dipoles in dilute systems is well known flij ] . Since 
we are interested in aqueous solutions, we have chosen 
as a fit parameter the molecular dipole moment of wa- 
ter to be po = 4.86 Debye (instead of the physical value 
pa = 1.85). This allows us to obtain e° = 80 for e = eo 
(vacuum permittivity) and Cd — 55 M. 

When the dipolar effects are strong (see below) there 
is a crowding of dipoles and ions between the plates, and 
their densities can reach values higher than close packing. 
To avoid this problem, we can generalize our theory to 
take into account the finite molecular size [l5| . Assuming 
that the 1:1 ions and dipoles are constrained on a lattice 
of spacing a (roughly equal to their molecular size), and 
imposing the condition that each site of the lattice is oc- 
cupied by only one of the three species (incompressibility 
condition), the free energy becomes 



pF= ^/d 3 r[W(r)] 5 
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where Cd + 2c s = a~ 3 . Minimizing the above free en- 
ergy, the Modified Dipolar Poisson-Boltzmann (MDPB) 
equation is obtained 
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The presence of the denominator V in (|12p leads to sat- 
uration of the local ionic and dipolar densities, which is 
quite important close to charged boundaries. Without 
the dipolar effect po — 0, the MDPB equation reduces 
to the modified PB equation which also displays an ionic 
saturation effect because of solvent entropy [la ]. 

A large deviation of the DPB treatment from the stan- 
dard PB one may occur in the strong E field regime. Such 
a case is presented now by solving numerically eqs. ([6])- 
J7|) for a system composed of two planar surfaces located 
at z — ±d/2, with opposite surface charge densities =Ftr 
and with small amounts of 1:1 salt to avoid strong screen- 
ing effects. In this anti-symmetric system the potential 
at the mid-plane vanishes, while the electric field there is 
non zero. The DPB pressure, in turn, deviates substan- 
tially from its corresponding PB value due to the coupling 
between the dipole density and the non-zero electric field. 
This is in contrast with a symmetric planar system where 
the electric field vanishes at the mid-plane. 

Had we modeled the water solvent as dipoles in vac- 
uum (e = e ), the dipole density in the mid region (see 
Fig. 1) would have reached unphysical values above the 
close packing ones, because nothing in our model pre- 
vents over-crowding. In order to avoid this artifact we 
use a background of low dielectric solvent {e.g., e r = 4.5 
for ether) and treat explicitly the strong water dielectric 
response by the dipolar term in the DPB equation ([5]) . In 
this fashion the water bulk density is lower than its close 
packing value, yielding a dipole profile density which is 
higher than the bulk value but below the close packing 
one. Note that all other mixture enthalpic and entropic 
terms are not considered at present 

In Fig. 1 we present the DPB profiles for the dipole 
density and local dielectric constant between two charged 
plates with separation of d — 20 A. The figure shows 
a strong accumulation of dipoles between the charged 
plates leading to high effective dielectric constant. The 
profile of the dipole density (dashed line) is rescaled by 
its bulk value. It can be seen that in the surface vicinity 
(up to about 2A) , the density rises to above four times its 
bulk value due to the strong attraction with the charged 
surface. In the mid-region the density saturates at about 
1.4 times its bulk value. The corresponding local effec- 
tive dielectric constant (solid line in the figure) can be 
calculated from @. The profile resembles that of the 
dipole density. In rescaled units, it saturates at a value 
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FIG. 1: The DPB rescaled dielectric constant e eff (2)/eb u ik and 
the dipole density Cd(z)/cd profiles between two oppositely 
charged plates at separation d = 20 A. The surface charge 
density is a — =pe/50A 2 . The reservoir contains 1:1 salt of 
concentration c s = 10 -5 M and dipoles of density Cd = 10 M. 
The dielectric constant is rescaled with respect to its bulk 
value ebuik = 18.2. The profiles have a strong variation in the 
vicinity of the plates (up to 2 A) and then saturate to a value 
that is somewhat higher than their bulk values. 



of about 1.2 in the mid-region and reaches about 2.3 at 
the surfaces. 

Compared to a PB theory with the same bulk dielec- 
tric constant (fT0|) , which is taken as constant throughout 
the system, the DPB demonstrates strong deviations, not 
only in the surface proximity but also in its saturated 
mid-range value (for strong enough a and/or small d). 

The ionic concentration is much less affected by the 
presence of the dipoles. We have computed the ion den- 
sities as a function of the distance to the surface. Because 
of the different boundary condition, ([7|), the ionic density 
is strongly suppressed at the surface with respect to PB 
(to about half of its original value). However, it comes 
back to its PB value at distances as close as 0.5 A from 
the surface. 

In Fig. 2 we plot the relative osmotic pressure differ- 
ence (Hdbp — npB)/npB as a function of the surface 
separation d. The pressure is a global quantity, and is 
sensitive to the strength of the electric field throughout 
the system rather than to its value on the surface. As 
a result, Hdpb deviates strongly from Hpb for small d, 
while Hdpb ~ Hpb at larger separation. 

We have presented an analytical modification of the 
PB equation by including the dipole degrees of freedom. 
We calculated the correction to the potential, electric 
field and densities for a system of two oppositely charged 
plates (Fig. 1). The results are compared with those of 
the usual PB equation with an effective dielectric con- 
stant. We find that when the electric field is strong 
(poE w ksT), there are strong deviations from the PB 
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FIG. 2: The DPB calculated pressure IT as function of the 
inter-plate separation 10 < d < 140 A for two oppositely 
charged plates. All other system parameters are as in Fig. 1. 
We compare the DPB with the usual PB models by plotting 
their relative difference (IIdbp — IIpB)/ripB. In the inset the 
DPB pressure is plotted in units of 10~ 2 fc.gT/A 3 . 

model. The spatial dependence of the dielectric constant 
signals an ordering of the dipoles at the surfaces. This 
spatial dependence is also a signature of non-linearity in 
the dielectric response. The inter-plate pressure is sensi- 



tive to the value of the electric field at the mid-plane and 
can deviate considerably from the PB results for small 
enough separation and/or large surface charges (Fig. 2). 

The formalism presented here is general but was ap- 
plied at a mean-field level. The PB equation is analyt- 
ically modified by the dipole degrees of freedom. As a 
mean-field approximation it also has its usual limitations. 
First, it lacks correlation effects. In addition, it does not 
treat correctly the finite size of the ions and dipoles and 
the densities of both can reach unphysical high values in 
high E fields. The latter limitation can be remedied by 
including the hard core of ions and dipoles, (fT2")l , and will 
be published elsewhere (l2T |. 

The results presented in this paper can be verified ex- 
perimentally by using, for example, the Surface Force 
Balance (SFB) apparatus. Recent SFB experiments fist ] 
have been performed on asymmetrically charged surfaces. 
The range of inter-surface separations that we used can 
be explored using the SFB technique. What is needed, 
however, are careful studies of mixtures of different di- 
electric solvents in order to extract the dipole contribu- 
tion to the osmotic pressure. This systematic set of ex- 
periments may shed light on the short-range hydrophobic 
effect and hydration forces. 
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